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We examine ways to write the Choptuik critical solution as the evolution of scale invariant vari- 
ables. It is shown that a system of scale invariant variables proposed by one of the authors does 
not evolve periodically in the Choptuik critical solution. We find a different system, based on max- 
imal slicing. This system does evolve periodically, and may generalize to the case of axisymmetry 
or of no symmetry at all. 
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I. INTRODUCTION 

Scaling behavior, as first found by Choptuik H, o ccurs at and near the threshold of black hole formation in the 
gravitational collapse of many types of matter ffllC|] . For a one parameter family of initial data slightly above the 
threshold, the mass of the black hole scales like {p — p*y . Here p is the parameter, p* is its critical value and 7 is 
a constant that depends on the type of matter, but not on the family of data. For initial data slightly below the 
threshold, the maximum curvature scales like {p* — p) '^'^ where 7 is the same constant as in the black hole mass 
scaling law. For some types of matter, the critical solution {p — p*) has periodic self-similarity: after a certain 
amount of time, the metric and matter variables repeat themselves with the scale of space shrunk. For other types of 
matter, the critical solution has exact self-similarity: only the overall spatial scale changes with time. 

Most of the work done on scaling has assumed spherical symmetry. However, Abrahams and Evans have found 
scaling in the collapse of axisymmetric gravity waves ^ . In addition, perturbative work of Gundlach on fluid collapse 
assumes no symmetries at all. Thus (for certain types of matter) scaling seems to be a generic property of critical 
gravitational collapse. 

The scaling of black hole mass (and of maximum curvature) has been explained p3|-p^ subject to the following 
assumptions: (i) the critical solution is periodically self-similar or exactly self-similar and (ii) the critical solution 
has exactly one unstable mode. What remains to be explained is why the critical solution is periodically self-similar 
(or exactly self-similar). It has been noted ||l|,|l3|,^ that the property of periodic self-similarity bears a striking 
resemblance to periodicity of a limit cycle of a dynamical system. Correspondingly, exact self-similarity of a critical 
solution resembles a limit point of a dynamical system. Furthermore, one can turn "resembles" to "is" by writing the 
Einstein field equations as a dynamical system of scale invariant variables. 

To write Einstein's equations as a dynamical system requires a choice of foliation (time slicing) for the spacetime. 
Some foliations will be compatible with the periodic self-similarity and some will not. More precisely, periodic self- 
similarity means that there is a diffeomorphism C, and a number A such that C* {gab) = e~'^'^ gat- A foliation is a one 
parameter family of spacelike (or null) hypersurfaces S(t). For a periodically self-similar spacetime, call a foliation 
compatible provided that for each to there is a ti such that C(I](io)) = ^(ti). In more physical terms a foliation is 
compatible if for each slice, there is a later slice that is identical except for an overall change in the scale of space. 
For a given periodically self-similar spacetime, there are many foliations that are compatible. Simply choose any 
initial hypersurface S(io)- Carry it forward by the diffeomorphism C to form the hypersurface S(ti). For to < t < ti, 
interpolate the hypersurfaces in any smooth way. Then extend that family to the whole spacetime using C. and C~^- 
There are also many foliations that are not compatible. Thus compatibility is a property of the entire foliation, or 
equivalently of the prescription for the lapse. Given any initial slice, there is some compatible foliation that contains 
that slice. A dynamical systems explanation of Choptuik scaling requires a compatible foliation. 

Note that compatibility is a condition on the lapse and is independent of the choice of shift. However, one may want 
to choose the evolution vector field so that the diffeomorphism ( has the form t ^ t + A with the spatial coordinates 
unchanged. This requires both a compatible foliation and a condition on the shift. 
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Most of the work done on Choptuik scaling has assumed spherical symmetry. Here the metric can be put in the 
form (I 

ds^ = - dt^ + dr^ + dn^ . (1) 

Here a and a are functions of r and t; dil^ is the unit two-sphere metric and r is the usual area coordinate. The 
surfaces of constant t are orthogonal to the surfaces of constant r. Usually the condition a = 1 at r = is imposed 
so t is proper time at the position of the central observer. This foliation is compatible. 
A spherically symmetric metric can also be put in the form 

ds"^ ^ ~gdu{gdu + 2 dr) + dVt^ . (2) 

Here r and dil^ are as in equation (1) and g and g are functions of r and u. The "time" coordinate u is constant along 
outgoing light rays and equal to proper time at the position of the central observer. This foliation is compatible. 

A metric or matter variable is scale invariant provided that it is unchanged when gab — > k gab for any positive 
constant k. The quantities a and a in equation (1) are scale invariant: under a scale transformation, a and a have the 
same form with only the overall scale of r and t changed. Similarly, the quantities g and g in equation (2) are scale 
invariant. In spherical symmetry, the metric is determined by the matter. The type of matter treated by Choptuik 
is a massless, minimally coupled scalar field (j). Initial data for this system consists of the field and its normal 
derivative at a moment of time. The matter variables used in |^ are (r / a) d(j)/ dr and (r / a) d(j)/ dt. These variables 
are scale invariant. For the same physical system, in the coordinates of equation (2) the usual matter variable is 
h = d{r(j))/dr, which is scale invariant. Scale invariant variables have also been found for other types of matter. 

While compatible foliations are known in the spherically symmetric case, Choptuik scaling is a generic phenomenon 
that does not seem to depend on symmetry. Therefore, one would like to find a compatible foliation for the general 
case. The foliations of equations (1) and (2) depend, in an essential way, on spherical symmetry. The constant t slices 
of equation (1) are defined as orthogonal to the spherical area coordinate r. The constant u slices of equation (2) are 
the null cones of the observer left invariant by the spherical symmetry. Thus, there does not seem to be any natural 
way to generalize these foliations to the case of axisymmetry or of no symmetry. 

Thus, what should be done is to use a slicing condition that makes no reference to symmetry, find a set of scale 
invariant variables associated with that foliation, and check that the foliation is compatible. In practice, one can check 
compatibility on a case by case basis. That is, given a particular critical solution, say as the output of a numerical 
code, one can check numerically whether a particular foliation is compatible. In this paper, we will explicitly check 
compatibility only with the Choptuik critical solution: that is, the case of a spherically symmetric, minimally coupled 
scalar field. 

In reference Q a slicing condition was proposed that makes no reference to symmetry, and a set of scale invariant 
metric and matter variables associated with this foliation was found. However, it was not known whether this foliation 
was compatible. This paper provides an answer to that question. In section 2 we present the numerical methods 
used in this study and consider the foliation of reference We show that this foliation is not compatible with the 
Choptuik critical solution. In section 3 wc find a different system of scale invariant variables based on maximal slicing, 
and show that maximal slicing is compatible with the Choptuik critical solution. We then consider the question of 
whether maximal slicing is compatible in the case of axisymmetry or no symmetry. 



II. NUMERICAL METHODS 



Typically, in numerical relativity, one thinks of the slicing condition as a part of the numerical code that evolves the 
metric variables to find the spacetime. Thus, the prospect of analyzing the Choptuik critical solution using several 
different slicing conditions seems daunting. A numerical treatment of the Choptuik critical solution requires great 
accuracy and stability and a large range of spatial scales. It would be difficult to write a code that does these things 
well in an arbitrary slicing. Fortunately, this is not necessary. Given the data from a code that uses one slicing 
condition, one can produce numerically the foliations of other slicing conditions. 

We begin with numerical data for the metric and scalar field for the Choptuik critical solution in the coordinates 
of equation (2), as found in reference (using methods based on the work of Christodoulou and of Goldwirth 
and Piran [|l9| ). A spherically symmetric slice in these coordinates is given by u = u(r). To produce a slice in these 
coordinates, we must express the slicing condition as a differential equation for u{r) and then numerically solve that 
equation. For the t = const, slices of equation (1), V at is orthogonal to V ar. Therefore these slices are given by 
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To solve this equation numerically, for each u we begin at r = and integrate outward using a Runge-Kutta method. 
For each (r, u) value we find the value of g by interpolation using the values that we have from the output of the 
numerical code of reference . 

Now suppose that we wish to find the maximal slices. For a surface u — u{r), the extrinsic curvature is 



2g du 



2 du 



^ ^^Ag dr ^ r 



Where W is given by 



Thus the equation K — becomes 



1 / T,.^ „xxA dg gW^ cPu 

' dr {du/drf dr^ 
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For each u we begin at r = and produce each surface by a Runge-Kutta integration outward using equation (6). 
Again, for each (r, u) value, the quantities g, g and their derivatives are interpolated from the gridpoint values of these 
quantities. 

Now suppose that our slicing condition involves a dynamical evolution of the lapse. Given the lapse N, we evolve 
the slice, and given the conditions on the slice, we evolve the lapse. For a slice, u = u{r), the unit normal vector is 



ou 



- w 



1 



du/dr 



(7) 



Since we wish only to know how the slices evolve, we can without loss of generality choose the shift to vanish. The 
coordinates of the slice then evolve by 



NW 



(8a) 
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du/dr 



(86) 



Here an overdot denotes derivative with respect to the evolution vector field N n°- . 
For the system of reference H The lapse evolves by 



1 o 

N = - K 

O 



(9) 



The 
and 



Let (hab, Kab, Da) be respectively the intrinsic metric, extrinsic curvature and derivative operator of a slice, 
scale invariant metric variables associated with this foliation are hab = N~'^ hab, Kab = {Kab ~ [K/3] hab 

UJa = DaN. 

We choose as the initial slice, a t = const, slice generated using equation (3). The initial lapse is chosen to be 1. 
The slice and lapse are then evolved using equations (8a), (8b) and (9) to produce a foliation of the spacetime. 

We would now like to know whether this foliation is compatible. In a compatible foliation, the scale invariant 
variables repeat their behavior. Thus at corresponding slices, the scalar field as a function of r/^/N should be the 
same. To identify which slices are "corresponding," define (/)o to be the value of </) at r = 0. Then <j)o repeats itself 
and the slices that contain repeating values of (/jq are corresponding slices. 
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FIG. 1. The lapse is plotted at the initial time and two later times. Note the steep peak in the lapse at the latest of these 
times. 

Note from equation (9) that positive K leads to an increasing lapse and, in fact, can cause the lapse to blow up in 
finite time. Figure 1 shows the lapse on three different slices. Note the steep peak in the lapse at the latest of these 
times. Shortly after this time, the lapse blows up. In the sequence of slices from the begining to the blow up of the 
lapse, 00 does not even go through one cycle. Thus this foliation does not even exist long enough for a comparison of 
corresponding slices to take place. Therefore this foliation is not compatible. 

III. MAXIMAL SLICING 

A dynamical systems explanation of Choptuik scaling requires a compatible foliation and a set of scale invariant 
variables. For the Einstein-scalar equations, initial data on a slice E are (hah-, Kab,(t'i P) where hab is the intrinsic 
metric, Kab is the extrinsic curvature, (j) is the scalar field and P is its normal derivative, (recall that for ria the 
unit one-form normal to S we have hab — 9ab + naUb, Kab — ha'^^ c^b and P — n°'\I a<P)- If (Safc, 0) is a solution of 
the Einstein-scalar equations and c is any positive constant, then [c^gab, (j^) is also a solution. In this new solution, 
the unit normal one-form is now ctIq. It then follows that the initial data for this solution are [c^hab, cKab, 4>, P/^)- 
Since a maximal slice in gab is also a maximal slice in c^gab, it is natural to define the following equivalence relation 
among maximal initial data sets: a set (/lab, Kab, 4>i P) is equivalent to a set {hab, Kab, (f', P) provided that there is 
a positive constant c such that 



{hab, Kab, 4>, P) = {c^hab, cKab, (t>, P / c) 



(10) 



Our set of scale invariant variables is an equivalence class of initial data. Is this set complete? That is, does an 
equivalence class of initial data contain enough information to determine its evolution? The answer is yes: simply 
take a member of the equivalence class, evolve it, and on each slice take its equivalence class. It does not matter 
which member of the equivalence class we pick, since in each case the metric that we produce differs only by an overall 
constant, leading to the same maximal slices and the same equivalence classes. 

For the case of vacuum spacetimes, the scale invariant variables are simply equivalence classes of {hab, Kab)- For 
types of matter other than the scalar field, one must add scale invariant variables corresponding to initial data for 
the matter fields. 
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An equivalence class is a somewhat abstract variable. Is there a way to make these scale invariant variables more 
concrete? In some cases, the answer is yes. Suppose that K'^^Kab is bounded on a slice E (as will be the case if E is 

1/2 

asymptotically flat) and define A = supj^ (if ^i^ab) . Then under the transformation gab (? gat we have A ^ A/c. 
Now define 



(hab, Kab, 0, P) = (A2 hab, XKab, ^, P/X) 



(11) 



Then under gab — > gab the set (hab, Kab, 4>: P) remains unchanged. In fact, this set is isomorphic to equivalence 
classes of {hab, Kab, 4>i P)- Thus for concreteness, in the case where K°-^ Kab is bounded on each slice, we can use 
{hab, Kab, 0, P) as our set of scale invariant variables. 
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FIG. 2. The scalar field is plotted as a function of rescaled radius 
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FIG. 3. The rescaled time derivative of the scalar field is plotted as a function of rescaled radius 



We now consider whether maximal slicing is compatible with the Choptuik critical solution. The maximal slices are 
generated numerically as described in the previous section. Compatibility of the foliation means that on corresponding 
slices, 4> and P as functions of Ar are the same. Note that this is all that one needs to check, since in spherical symmetry 
the metric variables are determined by the matter variables. In figure 2, </) is plotted vs. Ar for three corresponding 
slices (those with (/)o halfway between zero and its minimum). These slices correspond to A = 22.475, 706.33 and 
22361. Note that the three different curves are essentially indistinguishable, so (f> is the same on corresponding slices. 
In figure 3, P is plotted for these slices. Once again, the three curves are essentially indistinguishable, so P is the 
same on all three slices. It then follows that maximal slicing is compatible with the Choptuik critical solution. 

This numerical study looks only at the case of the massless, minimally coupled, spherically symmetric scalar field. 
Nonetheless, we now argue that maximal slicing is always compatible in the case of spherical symmetry, and may be 
compatible in the case of axisymmetry or no symmetry. Let gab be a periodically self-similar metric with associated 
diffeomorphism (. Let E be a maximal slice. Since C*(5ab) — e^^'^ gab, it follows that C(^) is a maximal slice. 
Therefore there must be a compatible maximal slicing that includes E. If the metric is spherically symmetric, there 
is a unique spherically symmetric maximal slicing. Therefore, this maximal slicing must be compatible. In the case 
of axisymmetry or no symmetry, things are more complicated because maximal slicing is no longer unique. The 
condition on the lapse TV that preserves maximal slicing is 



DaD^N + N Rabu^'n'' + R - ^^^i? 







(12) 



Solutions of this equation are determined by boundary conditions. For gravitational collapse, the usual condition is 
that iV ^ 1 at infinity. Perhaps this condition leads to a compatible maximal slicing. The work of reference Q on 
scaling in axisymmetric vacuum collapse uses maximal slicing with this condition on the lapse and finds that this 
foliation is compatible. It therefore seems likely that maximal slicing with iV — > 1 at infinity in general leads to a 
compatible slicing. In any case, there is some maximal slicing that is compatible, even in the absence of symmetry. In 
this foliation, with our scale invariant variables, critical gravitational collapse is simply a limit cycle of a dynamical 
system. 
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